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A thrust stand calibration technique for use in testing repetitively-pulsed electric thrusters for in-space 
propulsion has been developed and tested using a modified hanging pendulum thrust stand. In the imple- 
mentation of this technique, current pulses are applied to a solenoidal coil to produce a pulsed magnetic field 
that acts against the magnetic field produced by a permanent magnet mounted to the thrust stand pendulum 
arm. The force on the magnet is applied in this non-contact manner, with the entire pulsed force transferred 
to the pendulum arm through a piezoelectric force transducer to provide a time-accurate force measurement. 
Modeling of the pendulum arm dynamics reveals that after an initial transient in thrust stand motion the quasi- 
steady average deflection of the thrust stand arm away from the unforced or ‘zero’ position can be related to 
the average applied force through a simple linear Hooke’s law relationship. Modeling demonstrates that this 
technique is universally applicable except when the pulsing period is increased to the point where it approaches 
the period of natural thrust stand motion. Calibration data were obtained using a modified hanging pendulum 
thrust stand previously used for steady-state thrust measurements. Data were obtained for varying impulse bit 
at constant pulse frequency and for varying pulse frequency. The two data sets exhibit excellent quantitative 
agreement with each other as the constant relating average deflection and average thrust match within the 
errors on the linear regression curve fit of the data. Quantitatively, the error on the calibration coefficient is 
roughly 1% of the coefficient value. 


I. Introduction 

T HE evaluation of the performance of any propulsion system requires the accurate measurement of thrust. While 
chemical rocket thrust is typically measured using a load cell, 1 the low thrust levels associated with electric propul- 
sion (EP) systems for in-space propulsion necessitate the use of much more sensitive measurement techniques. The 
problem is compounded when a pulsed thruster is tested, as it subjects the thrust measurement device to an unsteady 
force. Work by Cubbin et al., 2 Haag, 3 and Markusic et a l. 4 present pulsed thrust measurement results obtained using 
torsional thrust stands where the pulsed thruster applied an impulse bit (impulse per pulse) to a swinging-gate thrust 
stand arm. In all three works the thrust balance response was calibrated by imparting an impulsive force (on the time 
scale of thrust stand motion) to the movable thrust stand arm using a swinging mass. The impulsive force imparted 
was found in Ref. [3] by calculating the momentum transfer during the impact from the measured change in velocity of 
the swinging mass and in Refs. [2,4] by directly measuring the imparted force using a piezoelectric force transducer. 
Knowing the thrust balance response to measured impulsive forces permitted the determination of a pulsed thruster’s 
unknown impulse bit by measuring the response of the thrust balance. Typically pulsed thrust measurements have 
been performed in single-pulse mode (single-shot) only, though Haag also performed repetition-rate thrust measure- 
ments by pulsing a LES-8/9 pulsed plasma thruster (PPT) on a timescale much shorter than the response time of the 
thrust balance arm, measuring the average displacement of the arm during pulsing, and equating that displacement to 
a ’steady-state’ average thrust level. 3 

In this paper, we describe a non-contact pulsed thrust calibration technique that is capable of electromagnetically 
applying a series of impulsive pulses to a thrust balance. The technique is capable of applying a single-impulse 
force for single-shot thrust measurement calibrations or providing multiple pulses to calibrate for a thruster operating 
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in a repetition-rate mode. While the technique of performing data acquisition and analysis on a single-shot test is 
well developed, the problem of operating and performing measurements on a thruster in repetition-rate mode has 
not been examined in great detail, with only a few researchers measuring the average thrust stand deflection for 
repetitively-pulsed thrust measurements. The technique is tested using the Variable Amplitude Hanging Pendulum 
with Extended Range (VAHPER) thrust stand, 5 demonstrating pulsed thrust measurement capabilities on this modified 
hanging-pendulum thrust stand that has only been used previously to measure the thrust levels of steady-state electric 
thrusters. 6 ' 7 The non-contact method is more easily controllable than the swinging impactor method used in previous 
works, applying impulsive forces to the stand in the same manner as the thruster. We also demonstrate in this paper 
that the measurement of the average thrust stand deflection during repetitively -pulsed calibration or operation is almost 
universally appropriate for the measurement of the average thrust and impulse bit, which for repetition-rate pulsing 
will obviate the need in the future to measure and anaylze the time-resolved motion of the thmst stand as is presently 
used when performing single pulse thrust stand measurements. 2,4 

The outline for the rest of this paper is as follows. In Section II we present an analytical model for thrust stand 
response to an arbitrary number of pulses. Section III contains a description of the test hardware used in the present 
effort. Data demonstrating the sensitivity and impulsive thrust measurement capabilities of the VAHPER thrust stand 
are presented in Sect. IV. 


II. Analytical Thrust Stand Response Model 

The response of most thrust stands, whether they are of the hanging or inverted pendulum or the swinging-gate 
variety, can be modeled as a torsional spring-mass-damper system. If the angular deflections are small, the deflection 
of the thrust stand arm can be linearized. For a repetitive, impulsive forcing function the governing equation can be 
written as 

*-?) (1) 

where x(t) represents thrust stand displacement as a function of time, ui n is the natural frequency of the thrust stand, 
and ( is the damping ratio. On the right side of the equation a repetitive impulsive forcing function with the magnitude 
/bit / Weft is applied to the system to represent the pulsing thruster or calibration unit, where / bit is the magnitude 
of the impulse and m e s is the effective mass of the spring-mass-damper thrust stand system. The impulsive load is 
applied as a delta function, 5 with a pulse width much less than the thrust stand natural period. This load is applied for 
a total of k pulses at a pulse frequency/. 

Solutions to the governing equation can be found for three separate cases, £ < 1, £ = 1, and £ > 1, corresponding 
to the thrust stand being underdamped, critically damped, and overdamped, respectively. A short derivation of the 
solutions for each case for initial conditions x(t) = x(t) =0 obtained using Laplace transforms can be found in the 
Appendix, while the more general solutions for an arbitrary set of initial conditions is given below as 
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where H is the Heaviside function. 

As an illustrative example of the response, the solutions given in Eqs. (2), (3), and (4) with initial conditions 
a;(0)=0.5 mm and ±(0)=2.5 mm/s for a single pulse (fc= 1) are given in Fig. la, showing the typical damped sinusoidal 
behavior of the system associated with an underdamped, critically damped, and overdamped system. Unless otherwise 
noted, all modeled solutions are for /bit=10 mN-s, m e ft-4.1 kg, and w, ,=1.57 rad/s. The time response of the system 
for an arbitrary number of pulses is presented in Fig. lb for all three damping cases, again with illustrative initial 
conditions ai(0)=0.5 mm and x(0)=2.5 mm/s. We observe that after a brief transient period, the system begins to 
oscillates about a given average displacement value. The transient response for the stand is different for each damping 
case, but the long-term average displacement value about which the system oscillates is essentially unchanged in all 
three cases with nearly equivalent time-histories after the first few seconds. We also calculated the behavior for varying 
initial conditions and found that the long-term behavior of the system does not change for any initial conditions. 
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Figure 1. Response of the system described by Eq. (1) subject to the initial conditions x(0)=0.5 mm and 4(0) =2.5 mm/s for a) a single pulse 
and b) multiple pulses at a pulse frequency of 1 Hz. 


In the simulated response data we generated to demonstrate the efficacy of the pulsed calibration concept, the 
average applied force to the system described by Eq. (1) is equal to the impulse bit times the pulse frequency ( / bit • /)■ 
If the relationship between this average applied force and the long-term average displacement value is linear, then a 
graph of the two should reveal the linear coefficient k e g, which during testing would represent the calibration constant 
relating a measured average thrust stand deflection to the unknown impulse bit produced by a thruster. This spring 
constant k e g is related to the natural frequency and effective mass of the system as 


2 ^-eff 

■ 

m e s 


(5) 


The deflection value is the difference between the zero position and the average value of the deflection once the 
response asymptotes to the long-term behavior observed in Fig. lb. The average force can be adjusted by either 
varying /bit or changing /. In Fig. 2a, /bit is varied at a constant pulse frequency of 1 Hz to apply a range of forces 
to the system. The graph exhibits a linear relationship between the applied forces and the average displacements 
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for all damping cases (only one case is shown for clarity since all three lie on top of each other). In Fig. 2b, the 
average force is changed by varying the pulse frequency at a constant impulse bit of 10 mN-s. We observe that 
the relationship between the applied force and the average displacement again follows a linear relationship, but only 
above the level where the pulse frequency at least twice as great as the system’s natural frequency of 0.25 Hz, with 
the discrepancy being most pronounced in the underdamped case. Physically, as the pulse frequency is lowered to this 
level the oscillator completes more than half a cycle between pulses and crosses zero, yielding a deflection response 
that is chaotic, never reaching the long-term repetitive pattern observed in Fig. lb. The difference between the average 
value of the long-term deflection and zero was selected as the metric, as opposed to using the difference between the 
maximum or minimum long-term deflection value and zero, because the slope of the curves given in Fig. 2 using the 
average deflection metric best matches the value of fc e ff used in Eqs. (2), (3), and (4) to generate the deflection curves. 



Figure 2. Average force applied to the system as a function of the average long-term deflection for a) varying impulse bit at a pulse frequency 
of 1 Hz and b) varying pulse frequency at an impulse bit of 10 mN-s. 


III. Experimental Apparatus 

A repetition-rate pulsed calibration was performed using the VAHPER thmst stand shown in Fig. 3a and previously 
described in Ref. [5]. The thrust stand is a variation on the hanging pendulum design, employing an adjustable 
mechanical linkage system to convert horizontal deflection of the pendulum arm to amplified vertical deflection of 
a secondary beam. The vertical deflection is measured using an optical linear gap displacement transducer (LGDT), 
while a passive eddy-current damping system is used to attenuate oscillatory motion of the balance mechanism. For 
the present work, a simulated pulsed inductive thruster having an approximate mass between 20-30 kg was mounted 
on the end of the pendulum arm to provide a realistic response of the thrust stand. 

An impulsive force was imparted to the thrust stand pendulum arm using the apparatus shown photographically in 
Fig. 3b. In this apparatus, a power supply is connected in series with the solenoid and switched using a high-power 
MOSFET to repetitively apply current through the solenoid coil. A square-wave input pulse to the MOSFET yields 
an approximately square-wave current in the coil, which is fixed in position and produces a magnetic field that acts to 
repel the magnet. The magnet transmits an impulsive (on the time scale of thrust stand motion) force to the pendulum 
arm through a piezoelectric force transducer that records the applied time-varying force and possesses a resolution of 
501.6 mV/N. The force applied to the magnet, and consequently to the thrust stand pendulum arm, is transmitted in 
a non-contact manner through the magnetic field produced by the current flowing in the solenoid coil. The recorded 
force (integrated with respect to time to yield the impulse bit) and the measured deflection of the balance mechanism 
permit the experimental generation of calibration curves similar to those presented in Fig. 2. 

IV. Calibration Data 

Thrust stand calibration data for repetitive pulsing are presented in this section demonstrating the applicability and 
usefulness of the analysis method outlined in Sect. II. Raw low-speed data showing the time-history of thrust stand 
deflection (given in V and acquired in the linear range of the position sensor) during a repetition-rate pulsed calibration 
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Figure 3. a) Schematic of the VAHPER thrust stand, b) Photograph of the pulsed calibration apparatus mounted on the thrust stand. 


trial are presented in Fig. 4. These data were acquired for a repetition rate of 10 Hz with force applied for 50 ms, with 
active pulsing over 30 s to reach a steady-state value. We observe that that when pulsing is initiated at t=0 the thrust 
stand deflects and oscillates as expected for a spring-mass-damper system until reaching a near-steady value beyond 
roughly 10 s. Once the forcing is discontinued, the same damped sinusoidal trend is observed about the zero deflection 
point. The difference between the zero position and the quasi-steady long term deflection during pulsing is taken as 
the average deflection for a measurement. 

Corresponding to the low-speed data in Fig. 4 are high-speed measurements in Fig. 5a of a single pulsed force 
applied to deflect the thrust stand arm. Temporally, this pulse occurs near the end of the pulsing time (near f=30 s in 
Fig. 4). This force is integrated with respect to time to yield the impulse applied during the pulse as presented in Fig. 
5b, with the maximum value at approximately 50 ms taken as the impulse bit / bit- We observe a 333 Hz oscillation in 
the force measurement, which are real oscillations in the applied force caused by current oscillations in the inductive 
solenoid coil owing to impedance mismatches in the electrical system. After integration the effect of this oscillation 
on the impulse is approximately zero. The average force applied to the thrust stand during a pulse train is simply the 
product of /bit and the pulse frequency. 

Calibration data are presented in Fig. 6a for a pulse frequency held constant at 10 Hz while the current through the 
solenoid coil, and consequently the applied force per pulse, is varied from trial to trial. In Fig. 6b the pulse frequency 
is varied in each trial from 1 to 10 Hz. The pulse width is 50 ms for all the data presented in Fig. 6. We observe that 
the average thrust as a function of average deflection is quite linear for both variable applied force and variable pulse 
frequency. More importantly, the linear calibration constant k e g in both cases is the same, within the calculated error 
on the regression fit, demonstrating the calibration technique described in Sect. II and showing the ability of the test 
apparatus to apply a wide range range of average forces to the thrust stand. The calibration constant errors for these 
data are roughly 1% of the total value of A: e ft ■ Data was acquired in Fig. 6b only down to 1 Hz since below this level 
the stand never reached a steady-state deflection following the linear trend but instead exhibited chaotic behavior, as 
was illustrated by the low-frequency modeling results in Fig. 2b. 
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Figure 4. Low-speed data showing the time-history of thrust stand deflection (given in V) during repetition-rate pulsing at a rate of 10 Hz 
over a 30 s timespan. 
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Figure 5. High-speed data showing a) the measured force and b) the impulse bit for a single 50 ms pulse driving the response of the thrust 
stand as shown in Fig. 4. 


V. Conclusions 

We have developed a non-contact thrust stand calibration technique for use in testing pulsed electric thrusters for 
in-space propulsion. The force is applied in a non-contact manner to a pendulum-based thrust stand by pulsing current 
through a solenoid coil that produces a pulsed magnetic field pushing against a permanent magnet mounted to the 
thrust stand pendulum arm. The force on the magnet is transferred to the pendulum arm through a piezoelectric force 
transducer, providing a real-time measurement of the force applied in a pulse. Unlike previous pulsed electric thruster 
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Figure 6. Measured average force applied to the system as a function of the measured average long-term deflection for a) varying impulse 
bit at a pulse frequency of 10 Hz and b) varying pulse frequency. 


work that focused primarily on single-shot pulsed measurements, this technique described in this paper has been 
developed for use on pulsed thrusters operating in a repetition-rate mode. Dynamic thrust stand modeling shows that 
after an initial transient in thrust balance deflection the average deflection away from the unforced or ‘zero’ position can 
be linearly related to the average pulsed force. This yields a simple linear Hooke’s law relationship between average 
displacement and the thrust where the spring constant k e g is the same as the one governing the dynamic motion of the 
thrust stand. The success of this method relieves the requirement to perform a more highly resolved temporal deflection 
measurement during repetition-rate pulsed thrust measurements in cases where the deflection during pulsing reaches 
a steady-state value. The method is not valid for pulse rates below a given frequency as the deflection never reaches a 
steady-state value. 

The technique was tested using a modified hanging pendulum thrust stand that had previously been used to obtain 
high-resolution steady-state thrust measurements on low-power Hall thrusters. Calibration data were acquired first by 
varying the impulse bit at constant frequency by varying the solenoid current, and consequently the magnetic field 
repelling the permanent. A second set of calibration data were acquired by varying the average thrust by changing 
pulse frequency at a relatively constant solenoid current. The linear calibration constant k e s obtained from both 
methods matched within the calculated error on the linear regression fit, and the overall error on the constant was 
roughly 1% of the coefficient value. 
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Appendix: Solution to the thrust stand deflection initial value problem 


The following is a derivation of the solution to the initial value problem of Eq. (1) subject to the initial conditions 
x(t) = x (t) = 0. The Laplace transform of the governing differential equation: 
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results in 
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where H is the Heaviside function. Note that these equations match those given in Eqs. (2)-(4) for x(t) = x(t) = 0. 
The general solution given in those equations can be obtained by following the same steps while keeping the general 
initial position and velocity terms in Eq. (A2). 
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